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Abstract 

We discuss chiral symmetry breaking in the intersecting brane model of Sakai and 
Sugimoto at weak coupling for a generic value of separation L between the flavour D8 
and anti-D8-branes. For any finite value of the radius R of the circle around which the 
colour Z)4-branes wrap, a non-local Nambu-Jona-Lasinio (NJL) type short-range interaction 
couples the flavour branes and anti-branes. We argue that chiral symmetry is broken in this 
model only above a certain critical value of the 4-dimensional 't Hooft coupling and confirm 
this through numerical calculations of solutions to the gap equation. We also numerically 
investigate chiral symmetry breaking in the limit R ^ oo keeping L fixed, but find that 
simple ways of implementing this limit do not lead to a consistent picture of chiral symmetry 
breaking in the non-compact version of the non-local NJL model. 
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1 Introduction 

The Nambu-Jona-Lasinio (NJL) model (T\ provides an example of dynamical chiral 
symmetry breaking C\fSB) and fermion mass generation in a simple effective field theory 
setting. In the original model, the fermions in the four-fermi interaction were taken to be 
the nucleons. Interest in the model has endured for two main reasons: (i) It appears to 
give a rather accurate description of chiral symmetry breaking and its consequences for 
low-energy hadron phenomenology; (ii) Appropriately replacing the original nucleons by 
coloured quarks, the model can be argued to describe all of the low-energy physics of QCD, 
including the anomaly term E [31 14| [1 . 

Recently, versions of the NJL model have emerged in a string theory setting [|6l |71, 
involving intersecting brane configurations. One such configuration is the model of Sakai 

'These works give an argument based on Wilsonian RG and the confinement property of QCD for the emergence of 
NJL model for quarks from the underlying microscopic dynamics, including the correct anomaly term with a coeflicient 
proportional to the number of colours. For a review of applications of this model to QCD phenomenology, see [5] 
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and Sugimoto (SS) [[8]|, which involves a system of intersecting D4, DS and anti-Z)8- 
branes. The SS model has been very successful in reproducing many of the qualitative 
features of non-abelian chiral symmetry breaking in QCD. In this model, the 'colour' 
Yang-Mills fields are provided by the massless open string fluctuations of a stack of A^, 
Z)4-branes, which are extended along the four space-time directions and in addition wrap 
a thermal circle of radius R. At scales much larger than string length, the theory on 
the D4-branes is (4 -i- l)-dimensional pure Yang-Mills with coupling = {2n)^gJs of 
length dimension. In the strong coupling limit, g^Nc » 2nR, this stack of D4-branes 
has a dual description in terms of a classical gravity theory [9] with the background 
geometry of a Euclidean black hole. Flavour degrees of freedom [QUI [TTl Wit are provided 
by the massless open string fluctuations between the colour branes and the 'flavour' DS 
and anti-D8-branes, which intersect the thermal circle at points separated by a distance 
L < nR. Various aspects of chiral symmetry breaking in this model have been discussed 

in [1 [II 121 [H [Bl [la [nl [III [IE I2Q1 |2B 121 12^ 

It was pointed out in ^ that the brane configuration of the SS model decouples the 
scales of chiral symmetry breaking and confinement §. The additional parameter in the SS 
model (as compared to QCD) which makes this possible is the ratio L/R. The authors of 
[[6l argued that in the limit R ^ oo with L kept fixed (non-compact SS model), the efl'ective 
low-energy description of the SS model at weak coupling § is given by a non-local version of 
the NIL model, which breaks chiral symmetry spontaneously at arbitrarily weak coupling. 
This result is surprising in view of our field theoretic intuition, which would suggest that 
chiral symmetry is not broken at weak coupling and that there is a transition to the broken 
phase at some critical value. One might suspect that this unexpected result is connected 
with the absence of a mass gap in the non-compact SS model, which results in a long-range 
four-fermi interaction. One way to test this hypothesis would be to work with a finite, but 
possibly very large, value of R, which corresponds to a confining theory with possibly a 
small, but non-zero, mass gap. From general arguments, one would expect chiral symmetry 
to be broken in this model at a length scale of the order of or smaller than the confinement 

similar observation was made in |32J in a different context. 
^There are several parameters of length dimension in the SS model, viz. R, L, Is and g^. As discussed in f6l, the weak 
coupling limit is defined by the hierarchy of scales g^Nc « Is « L « R. In this parameter region, stringy eff'ects may be 
neglected and, as we shall see, a controlled treatment of the interaction between left and right-handed flavours, mediated by 
Yang-Mills fields, can be given. The condition L « R makes it possible to have a chiral symmetry breaking (length) scale 
which is much smaller than the confinement scale. 



3 



scale Q However, in the corresponding effective non-local NJL model which describes the 
flavour brane-antibrane interactions, could be associated with length scales even larger 
than the confinement scale. This is because the NJL model does not incorporate confinement. 
Now, the point is that for solutions which have ;^'SB scale larger than the confinement scale, 
an efi'ective local NJL model should be adequate. But this latter model shows a critical 
coupling for;^fSB. This argument suggests that an effective non-local NJL model for the SS 
model would show^SB only beyond a critical value of the coupling. 

In this paper we analyze ;t'SB at weak coupling in the brane configuration of the SS model, 
with a finite value of R and a generic value of L. This theory has a mass gap, which gives the 
scale over which the four-fermi coupling of the flavour branes extends. We obtain the leading 
order (in gauge coupling) approximation to the effective fermion action, including the exact 
contribution due to the Kaluza- Klein modes. We expand on the plausibility argument given 
above that in the resulting non-local NJL model chiral symmetry is spontaneously broken 
only above a certain critical coupling. We then verify this by obtaining numerical solutions 
to the gap equation derived from the effective fermi theory. The plan of this paper is as 
follows. In the next section we first briefly review the argument of for the emergence of 
local NJL model from QCD and then extend it to a non-local model. We use this in Section 
3 to derive the non-local NJL model as the leading approximation to the coupling of the 
flavour branes in the weakly coupled SS model. In Section 4, we discuss ;^'SB in the non- 
local NJL model. We derive the gap equation in the large A^^. limit and present numerical 
solutions which show that chiral symmetry is spontaneously broken only above a certain 
critical coupling. A discussion of the non-compact case is given in Section 5. We end with a 
summary in Section 6. 

"*It is generally believed that in a confining gauge theory the length scale associated with is of the order of or smaller 
than the confinement length scale. QCD is an example where the;^'SB scale is of the order of the confinement scale. Recently, 
it has been argued in ll20ll2Tl that in the strongly coupled SS model the;^fSB length scale can be smaller than the confinement 
length scale, depending on the value of L. 
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2 NJL model from QCD 



The Yang-Mills action for U(Nc) gauge group (indices a, b) and Nf massless quark flavours 
(indices a,P) isS 

5o = -^ / + J d'xr{x)y>'(id, + fA"^{x)y{x). (1) 
This theory confines and develops a mass scale A, given by 3 

A~me"^, = -1^(11//, -2A^^). (2) 

It is generally believed that at energies below the confining scale, an effective NJL model 
for quarks captures the dynamics of the theory. There is no systematic way of integrating 
out the Yang-Mills degrees of freedom from QCD to get an eff'ective fermion action. A 
scenario outlining how one might think about doing this was presented in 01. The basic 
point is that integration of Yang-Mills degrees of freedom would lead to eff"ective multi-quark 
interactions. The range of these interactions must be short, of the order of 1/A, because of 
confinement, and so at energies below A a local approximation would be adequate. The 
NJL interaction between gauge-invariant quark bilinears is the leading term compatible with 
gauge symmetry and global symmetries of QCD. 



2.1 Extension to a non-local NJL model 

In QCD it is generally believed that the mass scale M associated with;t'SB coincides with 
the confinement scale A. Suppose, however, we can deform QCD in such a way that the 
two scales are separated by some new physics (as in the SS intersecting brane configuration 
discussed in the next section). In this case, for studying ;t'SB we need the efli'ective four-fermi 
theory at energies larger than the mass gap A. If A << M, the energies of quarks involved 
in the four-fermi interaction are much larger than A. Because of asymptotic freedom, for 
energies much larger than A we can present a more precise derivation of the efl'ective 

^We use the following notations and conventions. The space-time metric is mostly minus. Our Dirac matrices and 
their Weyl representation are as given in |33|, in particular, the equations (3.41) and (3.42). We have used the notation 

(u = 0, 1,2,3) to label the four space-time coordinates. Also, f are hermitian generators of UiNc) in the fundamental 
representation. In particular, we will need the identity {f)ij{f)ki - jSndjk- 

^This is true for Nf < 1 lNc/2, which is easily satisfied in the large Nc and fixed Nf limit that we will be interested in here. 
Also, the mass scale m that enters in this formula should be taken to be the scale at which the input coupling is measured. 
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interaction. The leading contribution to the effective interaction comes from a one-gluon 
exchange approximation, which can be calculated exactly. The result § is 



' efT 



-f/ 



d'xd'y ^^,{x-y)r''{x)r{y) (3) 



where 



7"^(x) = (ql{,x)a^fqL{x) + ql{x)(^fqR{x)^, (4) 



and 



— — , k^^l?-kl-ie. (5) 



74^ ^ik.x 

is the Feynman propagator for a massless scalar I^. Using the Fierz identities given in 
equations (3.77) and (3.80) of Jp3ll and retaining only interaction terms between left and 
right-handed Weyl components □, the effective action (|3]) becomes 

5o eff = ^4 J d'x d'y Ao(x - y)Vqt{x)(f^iym^iiy)qli.x)l (6) 

The bifocal fermion products in square brackets are singlets of the (global) colour U(Nc) 
group and transform as (Nf,Nf) under the flavour U{Nf) x U(Nf) group. 

The above discussion may be summarized as foffows. The effective four-fermi theory 
resuiting from integrating out the gfuon degrees of freedom is 



where 



5oeff = ij d'x(ql"(x)a^d,ql(x) + ql"(x)cr>'d,q''j,(x)^ 

+gl J d\d'yGoix-y)[qfix)q^^(y)][qf(y)qlix)l (7) 



ro(x) = r 7^ e'' -' Go(k), (8) 



^As usual, to do the calculation one needs to fix a gauge. The calculation done here and in the next section uses the 
Feynman gauge. 

^This rather unfamiliar way of writing the Feynman propagator is convenient since on making a Wick rotation to Euclidean 
signature and setting e to zero, k becomes just the magnitude of the Euclidean 4-momentum fc^. 

^The four-fermi terms involving Weyl components of a single handedness are not relevant to the discussion of chiral 
symmetry breaking vacuum. Hence, these terms are not taken into account here. 
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Gaik) satisfies: 

• For k « A, Go{k) ~ constant. This is ensured by the property of confinement of 
the action (O, which leads to the generation of a mass gap. The constant has length 
dimension two and we may take it to be 1 /A^. In a sense, this provides a definition of 
the mass scale A for us. 

• For A << ^ << M, Goik) ~ 1/P. This follows from the asymptotic freedom property 
of the action ([T]) and the added new physics which separates the ;^fSB scale from the 
mass gap of ©• 

A simple example of a function Go(k) which satisfies these two properties is 

^o(^) = TT^- (9) 
F + A^ 

A cut-off scale of order M on (|9l) is understood. A more complicated function could be 
devised to take into account the running of the coupling. In any case, the process of 
integrating out gluon degrees of freedom in a confining theory is expected to give rise to 
a far more complicated effective fermion action than in the model given by Qi, ^ and 
However, one might hope that the essential features for studying qualitative questions about 
;^fSB are present in this model. In the next section, we will see that a very similar non- 
local NIL interaction between the flavour branes emerges as the leading approximation to 
the weakly coupled SS model. 

3 NJL model from weakly coupled SS model 

At scales much smaller than the string length, the dynamics of the weakly coupled SS model 
is governed by the action [1^ 

S = J d^^J^ dx\Flj^(x,x^)f + J d^xql\x)d^(idf, + fA''^(x,-L/2)jql{x) 

+ J d'xql\x)(T^{id, + fA;(x,L/2)y^(x). (10) 



'"There are several possible corrections to this low energy action. For g^N « Is corrections from string modes are small 
and may be neglected. Corrections from the string winding modes around the thermal circle may be neglected for R » Is. 
We will assume this to be the case in the rest of this paper. The low energy effective action also has possible terms that couple 
the fermions to the transverse scalars. However, since the scalars come with a derivative, their effect may be neglected at low 
energies. 
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Only the space-time components A^{x, +L/2)cj of the (4+1 )-dimensional U (Nc) gauge field 
Am{x,x^) interact with the massless Weyl fermions q"^j^{x). Substituting the Kaluza-Klein 
expansion 

oo 

Al{x, /) = At\x) + YjfM^^y"'^" + At'\x)e-"'' (11) 

n=\ 

in this action, we get 

5=5o + Si + --- (12) 

So becomes identical to the action ([T]), after identifying the gauge potential A^(x) of the latter 
with the zero mode A'l^^\x) and setting 

gl = gj/2nR. (13) 

Also, it is now natural to identify the mass scale m in ^ with l/nR since the 4-dimensional 
description breaks down beyond this scale. S i is given by 

5 1 = ^ U J d'x(--\d,AfHx) - dAf'Hxf + yAf"\x)\') 

CO ^ 

+ J] j d'x(jr(x)AfHx) + JT(^)Af'^\x)\, (14) 

n=l ^ 

where we have used the notation 

r/(x) = {qlix)d^fqdx)e'"'^'^'' + qlix)<T^fqj,(x)e-'"''^''^. (15) 



The dots in (^12)) represent cubic and quartic interactions of the gauge fields. These will not 
be relevant to the leading order analysis in the weak coupling limit discussed below. 

We have already discussed the integration of the massless gluon degrees of freedom from 
the action So- Integrating out the massive Kaluza-Klein modes from the action (fT4l) is a much 
simpler task. To leading order in the gauge coupling, the effective four-fermi interaction due 



"in addition to the notations and conventions listed in Footnote 5, we use the following conventions. We use x'^ to label 
the coordinate along the circle which the ZJ4-branes wrap. We choose the mid-point between the locations of the D8 and 
anti-D8-branes on the circle, which are a distance L apart, as the origin in x''. The values = +L/2 are then the locations of 
the D8 and anti-D8-branes on the ckcle. 
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to the exchange of these modes is given by 

^ 1 eif = -^4 J ^-^^ - y)K;(.^)JT(y\ (i6) 



n=\ 

where 



is the Feynman propagator for a scalar of mass |. Using the Fierz identities (3.77) and (3.80) 
of [l33ll and, as before, retaining only interaction terms between left and right-handed Weyl 
components, the effective action (fT6l) becomes 

d'x cos(^)A„(.^ - y)yql"(x)cf^(y)][ql^y)qKx)]. (18) 

Now, using the identity 1.445.2 of [l34ll . 

cos ns n cosh a(n-s) 1 



Zcos ns 



, . 2a sinhna 2a^ 



we get 



5 1 eff = ^4 J A Gi(x - y)[ql"{x)cfj,(y)][qR(y)ql(x)], (19) 



where 



J (2;r)4 



k sinh K7r« 



A few comments are in order: 

• The interaction in (fT9l ), (|20l) is exact for all values of L and R, since we have summed 
over the exchange of all the Kaluza-Klein modes. In particular, in the limit i? — » oo, 
keeping L fixed, we get 

cm ^ ^e-''^. (21) 
k 

• For finite R, howsoever large, the second term in Gi(k) cancels the singularity in the 
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first term in the limit k « l/nR. This is consistent with our expectation that the 
range of the effective interaction (fT9l) should be of order the Kaluza-Klein radius R. In 
fact, for fermion momenta much smaller than 1 /R, one can approximate this non-local 
interaction by a local NIL term, with small derivative corrections. 

• For ^ >> 1 /L, the second term on the right hand side in (|20l) dominates, giving rise to 
a potentially problematic short-distance interaction with a 'wrong' sign. However, this 
term is cancelled by the large k contribution to the total effective action coming from 
the zero mode action, Q and Q. The net result is that for ^ » l/L, Gi(k) has the 
behaviour given in (|2T]) . 

3.1 Non-local NJL from SS model 

Combining (fT9l) with Q), we get the total effective fermion action 

+gl J d\d'yG(x-y)[ql"(x)q^^(y)][qf(y)ql(x)], (22) 

where 

r d'^k ., ~ 

G{x) = -—e''-'G{k\ (23) 

and G{k) = Go{k) + Gi(k). Although a precise derivation of Go{k) does not exist, it must 
satisfy certain conditions which we discussed in the previous section. As a result of these, 
G(k) must satisfy: 

• For k ~ A, G(k) ~ 1 /(k^ + A^). This follows from the hierarchy of scales A « l/nR < 
1 /L and the fact that the range of the non-local quartic fermion interaction in (l22l) is 
set by the mass gap dynamically generated in the 4-dimensional Yang-Mills action ©, 
which is of order the glueball mass ~ A. 

• For A << ^ << l/nR, G{k) ~ 1/P. This is because for these values of k, the 4- 
dimensional description continues to be valid and we are in the asymptotically free 
regime. 

• For k » l/nR, G(k) ~ nRe'^^/k. Here we have assumed L « nR (which is 
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automatically satisfied in the limit of large R, keeping L fixed). In this regime, the 
4-dimensional description is inadequate since the Kaluza-Klein states can be easily 
excited. We must now use the full 5-dimensional description. The exponential fall- 
off implies that the quartic fermion interaction in (l22l) has a short-distance cut-off as 
well. The effective scale here is of order L, the separation between flavour D8 and 
anti-Z)8-branes. 

A simple function that contains all the three scales and seems to capture all the essential 
features discussed above is 



1 +nRk 



This differs from ^ in that 1/L enters as a smooth ultraviolet scale, as opposed to the hard 
cut-off M that came with @. Moreover, it contains the additional scale R, which has to do 
with the underlying 5-dimensional origin of the model. We will use this simpler function in 
the analysis that follows. Also, throughout the following we will assume L « nR. 



4 ;^fSB in non-local NJL model 



As usual, we first introduce scalars to rewrite in the equivalent form 

Seff = /J d'x((i^{x)a^d,ql{x)) + q'^{x)(T^'d,ql{x) 

T"P\x,y)T"l^{x,y) 



+ T">'^(x,y)qf(y)qlix) 



glG{x-y) 



+r^x,y)ql"{x)4(y) 



(25) 



It is easy to verify the equivalence of this action to (|22|) by using the equations of motion of 
the scalars. 



T"^x,y) = gl G(x - y) qf(y)qUx). 



(26) 



The next step is to integrate out the fermions to get an effective action for the scalars. In the 
\arge-Nc limit, a classical treatment is adequate. Since we are only interested in the solution 
corresponding to the ground state, which is Poincare invariant and invariant under diagonal 
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(vector) flavour group, we may use the ansatz T"'^{x,y) = 5"^T{\x - y\). This simplifies 
calculation of the effective action. Making a Wick rotation to the Euclidean signature, we get 



VN,.Nf g 



Lr,.,E«t_ r^Wi.rat), (27) 



where V is the 4-volume and T{k) is related to T{x) by a Fourier transform, which we define 
by the relation ([8]) for any function. Also, G^(x) is the Euclidean version of G{x). Taking the 
Fourier transform of (|24l) . we get 

where k is the magnitude of the Euclidean four-momentum and Ji(k\x\) is a standard Bessel 
function. The ^-integral can be done using the identity 

l+nRk 1 + inRA n , ... , 

——- = ——— dse"''e-''+c.c., (29) 
k^ + 2iA Jo 

and the identity 6.623.2 of [l34ll . After some simplification, the final expression takes the 
form 

G'^ix) = ^ g{\x\A), (30) 

where 

= TP) — ^7TT77 + (cosLa +i?AsinLA) Ji(r) + (sin La - i?A cos La) I^ir). (31) 
In the above, we have introduced the dimensionless quantities 

i?A = TT^A, = ^A. (32) 

Moreover, 



^^'^ Jl. (*' + ^2)^^/^ r Jo -{s'^r^W' 



COS 5 



ds- 



r°° sin 5 1 TT sin 5 

^2(r) = ds =---(/i(r)-Li(r))- (33) 

(5^ + H)J/2 r 2r Jo (i'^ + r^yl^ 
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where Ki(r) and /i(r) are standard Bessel functions and Li(r) is a Struve function. It turns 
out that numerical calculations are done faster with the second form of 1 12(f). Note that in 
the region r « i?A, all the 5 dimensions must come into play. Indeed, in this region we get 



g(r) 



Ra/LI for r«LA 
i?A/r^ for La « r « Rf 



(34) 



For r » Ra, the system becomes effectively 4-dimensional. Here we find 



gir) 



l/r^ for 1 » r » i?/ 



(35) 



In Figure [U we have plotted the function g{r) as a function of r. We see that it approaches 
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50000 
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0.08 



Figure 1: g{r) as a function of r for the parameter values La 
Ra = 0.2. 



0.01 and 



a constant for r « La- For r > La, it decreases very rapidly and eventually for r > 1 (far 
beyond the region shown in the figure) it decays exponentially!^. 



'^This behaviour actually holds only for intermediate values of r which satisfy r'l^e^'' > Ra- If Ra is small, this inequality 
can allow rather large values of r. At much larger values of r, g{r) decays only as l/r^. This is presumably an artifact of 
the choice of the G(k) function we have made in (l24l i. For practical reasons, in the numerical calculations done in the next 
section, we have simply set g(r) to zero beyond a sufficiently large value of r. This is consistent with the expectation that the 
interactions should decay exponentially beyond the confinement scale, r = 1 . 
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4.1 Gap equation and order parameter of ;^fSB 

The equation of motion for T(x) following from the effective action (ITTl) is 



d'^x -^e-""-' = (36) 



This nonlinear equation for the order parameter is the analogue of the gap equation for the 
present case. T{x) = is the trivial solution which preserves chiral symmetry. However, this 
is not the solution which minimizes the effective action (|27l) . To see this multiply (l36l) 
on both sides by f*(k) and integrate over k. This gives 



Using this in (|271) . we get 



%J G^ix) J i2nYk^ + \f(kW ^ ^ 



sf, r d'k 



VNcN 



f 



^ r d'^k r \T{k)Y I ^ \T{k)\' \ 



(38) 



It is easy to see that the integrand on the right-hand side above is a decreasing function of 
\t{k)\lk and that it vanishes for T{k) = 0. It follows that T(k) = is not the solution which 
minimizes the effective action (1271 ). We also note that a potential divergence from the large k 
end gets cancelled between the two terms in the integrand and the net result of the integration 
over k is finite, provided f{k) is a decreasing function for large k. If such a solution exists, 
then the chiral symmetry is spontaneously broken. 

The order parameter of chiral symmetry breaking is the condensate 

m = ^{ql"(x)q20)}. (39) 

1\ r 



The field T{x) is related to it by UB, i-e. 

T{x) = 4n^A G^{x) (p{x), A = glNJAn^. (40) 



The gap equation (1361) can be rewritten in terms of (p{x) as 



= . (41) 

F + \T{k)\^ 



14 



We will look for solutions of this equation with ^(jc), and hence T(x), real. Since these are 
spherically symmetric functions of \x\, their Fourier transforms, (pik) and f(k), are also real 
functions of k. Furthermore, we see from (|4T]) that for large k, (p{k) « l/k. This is because 
f(k) must be a decreasing function for large k, for reasons explained above. Now, solving 
(|4T]) for f(k) as a function of 4>(k), we get 



f,ik) ^ 



1 + 



^1 - 4F02(/t) 



(42) 



This has real solutions only for 0(^) < 1 1 2k. For large fc, with 0(fc) << 1/fc, we get 
T+{k) » k and f_{k) « k. So, the desired solution f{k) must coincide with f_{k) for 
large fc. For small enough k, the solution f(k) must go to a constant (the mass gap) and so 
it must coincide with f+(k) for small k. The transition from one to the other occurs at some 
scale k = IX, where the two solutions coincide, i.e. f+{jx) = T^(ju). From (|42l) we see that /i 
satisfies the equation ^(ju) = l/lju. 

4.2 Solutions of the gap equation 

We parametrise the condensate as follows: 

m = ^<^(ui//). (43) 

Here / is the ;^SB scale. The normalization has been chosen to explicitly display the 
dimensions of the order parameter and to have i^(0) = 1. Now, using (|40l ). the Fourier 
transforms, ^(k) and T(k), can be written as 

1 r°° 

m = Icpo m\ f(p) ^- dy /Mpy) ip(y), 

P Jo 

1 r"' 

T{k) = AUh(t>ot{kl), m = - dyy^Mpy)g{Uy)(p(y), (44) 

P Jo 

where /a = ZA and Ji is a standard Bessel function. Using these, we can rewrite the gap 
equation (|4TI) as 

= 2 t^!2 2. ^ - ^ = -^^A- (45) 
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This is a nonlinear equation which we are unable to solve analytically. However, there are 
some general observations we can make. 

• Equation (|45l) cannot have a solution with the;^fSB scale / arbitrarily smaller than L. On 
physical grounds, we expect ip{\x\ll) to be substantially different from zero only in the 
region < |;c| < /, vanishing rapidly for \x\ » I. As a result, most of the contribution to 
the integral in the definition of t{p) in (|44l) comes from a region in which the argument 
of the function g{r) varies over the range [0, /a]. For I « L, g(r) is roughly a constant 
(= ff) over this range of r, as can be seen from (|3T]) . So, for / « L, t(p) » /(p). But 
this is not a solution to (|45|) . as can be easily checked. This argument works even better 
for large p, because then the contribution to the integral from y > I region is even more 
suppressed, beyond that due to a rapidly falling (p(y). But, for large p a solution to (l45l) 
must satisfy t{p) <x p^f(p). Thus, the gap equation (l45l) has no solutions for / << L. 

• In principle, ;(fSB solutions to the gap equation should exist for all l> L. This is because 
our non-local NIL model does not incorporate confinement. However, for solutions 
which have / > A~\ i.e. /a > 1, it should be possible to replace the non-local model 
by an effective local NIL model with A as the ultraviolet cut-off. This is because of 
the exponential decay of the four-fermi interaction g{r) for r > 1 . Since the local NIL 
model has no ;^'SB solutions below a critical coupling, we expect a critical coupling to 
show up for values of order l/>^> I'm the present non-local NIL model as well. One way 
this can happen is that as I a increases beyond La, the value of A which gives rise to this 
solution decreases until it hits a critical value at around Za ~ 1 . This argument cannot 
be made for the non-compact SS model considered in |l6l because of the absence of the 
mass scale A and the consequent absence of the exponential decay of g{r) for r > 1. 

Numerical calculations reported below bear out both the above expectations. 

4.3 Numerical solutions 

In the following, we will report on some solutions to the gap equation (|45|) obtained 
numerically using mathematica. This numerical work is based on the following strategy. 
From the expected physical properties of the order parameter (f>(x), we first make an ansatz 
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for it: 



The power cr and the constant c are adjustable parameters. With the scale / and the 
normalization constant cpo, there are altogether four adjustable parameters. Given (|46|) . the 
two sides of (l45l) can be computed and compared, and the difference can be minimized by 
varying these parameters. The numerical computations were done as follows. For the left- 
hand side of (|45|) . we need to calculate fip). This can be done once we choose some values 
for cr and c. After some experience with the calculations, it was not hard to make a good 
guess for the right values. For the right-hand side of (|45l) . we need t{p) for which one first 
needs to calculate the function g{r). This requires choosing values for i?A and La. In Figure 
[Uwe have shown an example of g(r) for Ra = 0.2 and La = 0.01. With g(r) at hand, one 
can now calculate t(p), after making a choice for /a. The right-hand-side of (|45] ). which we 
denote as ft(p), 

''^'^ -MP), (47) 



p^ + A^(f>lf(p) 



can then be computed. This requires making a choice for the parameters A and ^q, which 
are adjusted such that the deviation J dp(f(p) - ftip))^ is minimized. In principle, in the 
calculation of the deviation, the range of the integral over p should extend to infinity. In 
practice, we have found a value of about 10 to be good enough for the upper limit (for the 
values of the parameters La and we have used in our calculations), in the sense that the 
value of the integral remains essentially unchanged if the upper limit is increased beyond this 
value. The whole procedure was then repeated with slightly different values of cr, c, A and 
00 until the deviation was minimized for the chosen values of R^, La and /a. The value of A 
relevant to these values of the parameters was obtained from its relation to A given in (|45l) . 

In Figure [2] we have given two examples of the quality of solutions obtained in this way 
for two different sets of values of {i?A, La}, with Za >> La. The agreement between f{p) 
and f,{p) is excellent. In fact, the deviation J dp(f(p) - ftip))^ is less than a thousandth of a 
percent of the quantity J dp(f(p))^. Surprisingly, in all the solutions that we have obtained, 
the value of the parameter c which minimizes the deviation turns out to be exactly l/L. This 
may provide a hint for analytical solutions of the gap equation. 
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p 



Figure 2: The functions f(p) and ft{p) shown in these figures are the two 
sides of the gap equation. The two figures correspond to two different sets 
of values of {R^, L^}, as indicated. 

Note that both cases in Figure [21 have /a » La. As argued in the previous subsection, 
we do not expect any solution for Za << La. In fact, our numerical calculations indicate that 
there is no solution for 1^ < L^. This can be seen from the example given in Figure [3] where 




Figure 3: The two sides of the gap equation for /a = La/5. 

we have taken I a = La/5. There is no agreement, which is the best we have been able to do 
with the ansatz (|46l ). The seeming agreement at large p is misleading because the values of 
both f{p) and /,(/>) are so small that the figure cannot distinguish them from zero at the scale 
used. A better measure for the behaviour at large p is the ratio t(p)/f(p), which is expected 
to approach a constant at large p for Za << La. In the first of Figured we have plotted this 
ratio. It behaves as predicted at large p. The value of the constant also turns out to be close to 
the expected one, namely 7?a/L^. For comparison, in the second of Figure S] we have plotted 
the ratio p^f(p)/t(p) for the second example of Figure |2l which is expected to approach a 
constant at large p since this provides a solution to the gap equation. The figure verifies this, 
implying that in this case t{p) ~ p^f(p) at large p. 

In Figure |5] we have plotted Za as a function of A. Figures for two different sets of values 
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f(p) t(p) 




p 



Figure 4: The first figure shows t(p)/f(p) for Za < La and the second figure 
shows p^f(p)/t{p) for /a >> La as a function of p. 



of {Ra,La} have been given. Starting with Za = La, we see that at first A decreases with 



La-- 



300 



Figure 5: Za as a function of A for two different sets of values of {Ra, La). 



increasing Za, until it reaches a minimum at around Za = 1. Beyond this point, increasing 
Za seems to be accompanied by an unchanged or perhaps even an increasing A\}j. We have 
verified the behaviour shown in the two examples in Figure [5] for several other values of 
the set {Ra, La} and believe that this is the general behaviour. These data provide fairly 
convincing evidence for the existence of a critical value of A below which no solutions to the 
gap equation exist. 



5 The non-compact limit 

It is of interest to ask what happens in the non-compact limit, R ^ oo keeping L fixed. In our 
model, taking this limit is somewhat subtle, because the hierarchy of scales L « nR « A"^ 

'^This point is difficult to clarify with much accuracy beyond the range of values of Za shown in the figure. This is because 
calculations for such large values of /a require a greater precision in calculating g{r) and hence take much longer time. 
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must be maintained as the limit is taken. Furthemiore, even though the non-local NJL model 
does not incorporate confinement, to connect with the underlying weakly coupled gauge 
theory we may wish to impose on this model the relations between the confinement scale 
A~\ R and the coupling g^, given by ^ and ([13]) (with m ~ l/nR). Under the scaling 
R — > R(ri) = rjR, the following scaling properties can be deduced from these relations: 

Ra(7i) = {Ra)\ Uiri) = La^^^, (48) 

where La and i?A are the values for j] = I. The non-compact limit corresponds to taking 
?7 — > cx). Since Ra « 1, this means that both i?A('7) and La(?7) vanish in this limit. 

Let us denote by the value of A for I a = 1 . A^ is close to the minimum value of A 
(see Figure [5]) and so may be considered to be the value of the critical coupling. How does 
Ac change as a function of ?7? To find this, we have numerically calculated /If (77) for the set 
{7? A (?7), La (77)} for different values of 77. The numerical data have been plotted in Figure [6] 
as a function of 1/77. Calculations were done for two different sets of values of {Ra,La} to 
check dependence on initial conditions. Numerical data in the two graphs of Figure [6] show a 
similar pattern, indicating no dependence of the general behaviour on initial conditions. The 




0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1.0 



Figure 6: Ac as a function of I/77. The data in the two figures correspond 
to the scaling rule (1481) for two different initial values of the set {Ra,La}, 
{1/2, 1/100} for the first figure and {1/2, 1/30} for the second. The solid 
lines are drawn to indicate the region in which linear behaviour with I/77 is 
seen. The general behaviour is similar in the two cases. 

data show that at first, as 77 grows away from the initial value 77 = I, Ac decreases linearly with 
1/77, as one might expect from the relation between the 4-dimensional 't Hooft coupling and 
its 5-dimensional counterpart, A = As/lnR, which was used in deriving the scaling relations 
(ffSl ). However, at some point, A stops falling and seems to bottom out and start increasing or 
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perhaps go to a constant o This is not consistent with the relation between the couplings, 
so it seems that the non-compact limit cannot be reached by the one-parameter scaling given 
by®. 

An alternative way of taking the non-compact limit is as follows. The scaling rule 
(1481) used in the above analysis was derived from relations between A, R, A and A5 which 
follow from confinement in the underlying low-energy gauge-theory. Since the non-local 




Figure 7: Ac as a function of I/77. These data have been obtained for 
La{J]) = L/^/j], with La = 3/200 and keeping 7?a fixed as 77 is increased 
from 1 . The solid line is a power-law fit to the data. 



NJL model does not incorporate confinement, we may wish to relax these conditions among 
the parameters. However, we must still maintain the hierarchy of scales L « nR « A'^ . 
One simple way to do this is to keep R\ fixed (and << 1) as i? is scaled. Under the scaling 
R — > riR, then, we must have A ^ A/r]. This implies that La scales to La///. 

In Figure |7] we have shown data for Adr]) obtained as a function of 1 /r] by computing the 
critical coupling for the set {Ra, LA(ri)} for different values of 77. The data fit almost perfectly 
to a power law, Ac = 0.271339/77*'^'^'^^. Since this fit implies that TjAdrj) blows up in the 
limit 77 — > 00, the non-compact limit cannot be reached by this scaling either. 

Our tentative conclusion is that possible solutions in the non-compact version of 
the non-local NJL model cannot be obtained by taking any simple R ^ 00 limit of the ;^'SB 
solutions of the compact model. Clearly this issue deserves to be investigated further. 



'^More detailed calculations are needed to settle between these two possibilities. Calculations at higher values of 77 involve 
fine-tuning of various parameters of the solution to the gap equation and so they are harder to determine. For this reason we 
have restricted ourselves to a maximum of 77 = 5.5. It should be possible to go to larger values with some more effort. 
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6 Summary 



The study of ;fSB in QCD is made complicated by the fact that the scale at which chiral 
symmetry is broken is of the order of the confinement scale. If QCD could be deformed to 
enable "tuning" of the scale to be much smaller than the confinement scale, then one 
would have separated the complications of the dynamics of confinement from a study of ^SB, 
which could then be handled by perturbative methods. The intersecting brane configuration 
of Sakai and Sugimoto, which gives rise to a QCD-like theory at low energies, admits just 
such a possibility; it has an additional parameter, the flavour brane- anti-brane separation, 
which can be tuned. In the strong coupling limit, by tuning this parameter one can indeed 
raise the chiral symmetry restoration temperature above the deconfinement temperature 
[|20ll2T]| . However, as we have seen in the present work, for a consistent description of;t'SB 
in the weakly coupled SS model, it is essential to incorporate the physics of confinement. 
The interaction between flavour branes and anti-branes of the SS model is captured by a 
non-local NIL model. For any finite radius R of the circle which the colour D4-branes wrap, 
there is confinement and a mass gap in the low energy theory. The NIL model reflects this 
dynamically generated mass scale. A, in the length scale over which the non-local four-fermi 
interaction extends. This fact turns out to be crucial in getting consistent SB solutions. In 
the large Nc limit, the question of ^SB amounts to finding appropriate solutions to the non- 
linear gap equation. For solutions with;^SB length scale / much larger than the confinement 
scale A"\ it is reasonable to replace the non-local NIL model by the local NIL model. Hence 
these solutions must reveal the existence of a critical coupling, which is known to determine 
^SB in the local NIL model. In this paper we have numerically solved the non-linear gap 
equation and verified the existence of a critical coupling below which chiral symmetry is 
unbroken. Roughly speaking, only solutions with ;^SB scale greater than the brane-anti- 
brane separation L exist. The ;^SB scale I increases as the 't Hooft coupling is decreased, 
until a critical coupling is reached for / ~ A Solutions with / > A Mo not lead to any 
further decrease in the coupling. 

Our analysis is valid for any finite value of the radius R, which may be large. We have 
briefly addressed the question of what happens when 7? — > oo. Two diff'erent ways of taking 
this limit, each one obtained from a well-motivated one-parameter scaling of the parameters 
of the SS model, were discussed. We found from our numerical data that neither of them 



22 



leads to a sensible limit. The tentative conclusion is that simple ways of implementing this 
limit do not lead to a consistent picture of ;^SB in the non-compact version of the non-local 
NJL model. This seems to reinforce the critical role that the confinement scale plays in 
the compact model; the infrared cut-off provided by it enables the existence of consistent 
solutions to the gap equation. However, more work needs to be done to clarify this issue 
further. 

Finally, most of the calculations reported in this paper were done numerically because 
the gap equation is non-linear and we could not solve it analytically. It would, however, be 
useful to have some analytic handle on the calculations, especially in the parameter region 
near the critical coupling. This could be important for a better understanding of the non- 
compact limit. A possible hint in this respect is the fact that excellent numerical solutions 
were obtained using the ansatz (|46l) . with the constant c turning out to be almost exactly 
equal to Z/L in all cases. 

References 

[1] Y. Nambu and G. Jona-Lasinio, Dynamical Model Of Elementary Particles Based On 
An Analogy with Superconductivity. I, Phys. Rev. 122 (1961) 345. 

[2] A. Dhar and S. R. Wadia, Nambu- Jona-Lasinio Model: An Effective Lagrangian for 
Quantum Chromodynamics at Intermediate Length Scales, Phys. Rev. Lett. 52 (1984) 
959. 

[3] A. Dhar, R. Shankar and S. R. Wadia, Nambu- Jona-Lasinio-type effective Lagrangian: 
Anomalies and nonlinear Lagrangian of low-energy large-N QCD, Phys. Rev. D31 
(1985) 3256. 

[4] A. Dhar, Chiral symmetry, large-N limit and effective lagrangians in strong interactions. 
Invited lecture delivered at the University of Delhi Workshop in Physics (Jan. 1984), 
TIFR report TIFR/TH/84-25. Unpublished. 

[5] T. Hatsuda and T. Kunihiro, QCD phenomenology based on a chiral effective 
Lagrangian, Phys. Rept. 247 (1994) 221,jhep-ph/940B10, 

[6] E. Antonyan, J. A. Harvey, S. Jensen and D. Kutasov, NJL and QCD from string theory, 
hep-th/0604017| 

23 



[7] E. Antonyan, J. A. Harvey and D. Kutasov, Chiral symmetry breaking from intersecting 
D-branes, Nucl. Phys. B784 (2007) 1, hep-th/0608177, 

[8] T. Sakai and S. Sugimoto, Low energy hadron physics in holographic QCD, Prog. 
Theor. Phys. 113, 843 (2005), hep-th/04 12141] 

[9] E. Witten, Anti-de Sitter space, thermal phase transitions, and confinement in gauge 
theories, Adv. Theor. Math. Phys. 2, 505 (1998), ,hep-th/9803T3T ] 

[10] A. Karch and A. Katz, Adding flavour to AdS/CFT, JHEP 0206 (2002) 043, 
Lhep-th/0205236, 

[11] J. Babington, J. Erdmenger, N. J. Evans, Z. Guralnik and I. Kirsch, Chiral symmetry 
breaking andpions in non-supersymmetric gauge/gravity duals, Phys. Rev. D 69 (2004) 
066007, hep-th/03060T8l 

[12] M. Kruczenski, D. Mateos, R. C. Myers and D. J. Winters, Towards a holographic dual 
oflarge-N(c) QCD, JHEP 0405 (2004) 041, |hep-th/0311270| 

[13] T. Sakai and S. Sugimoto, More on a holographic dual of QCD, Prog. Theor. Phys. 114, 
1083 (2006), hep-th/0507073, 

[14] K. Nawa, H. Suganuma and T. Kojo, Baryons in holographic QCD, Phys. Rev. D75 
(2007) 086003, hep-th/06 12187) 

[15] K. Nawa, H. Suganuma and T. Kojo, Brane-induced Skyrmions: Baryons in 
holographic QCD, Prog. Theor. Phys. Suppl. 168 (2007) 231, |hep-th/0701007l 

[16] D. K. Hong, M. Rho, H. U. Yee and P. Yi, Chiral dynamics of baryons from string 
theory, hep-th/0701276. 

[17] H. Hata, T. Sakai and S. Sugimoto, Baryons from instantons in holographic QCD, 
rhep-th/0701280, 

[18] O. Bergman, G. Lifschytz and M. Lippert, Holographic nuclear physics, 
[arXiv:0708.0326) 

[19] D. Yamada, Sakai-Sugimoto model at high density, larXiv:0707.0101l 

[20] O. Aharony, J. Sonnenschein and S. Yankielowicz, A holographic model of 
deconfinement and chiral symmetry restoration. Annals Phys. 322 (2007) 1420, 
|hep-th/0604161. 



24 



[21] A. Parnachev and D. A. Sahakyan, Chiral phase transition from string theory, Phys. 



Rev. Lett. 97 (2006) 1 1 1601, hep-th/0604173 



[22] R. Casero, E. Kiritsis and A. Paredes, Chiral symmetry breaking as open string tachyon 
condensation, Nucl. Phys. B787 (2007) 98, hep-th/0702155. 

[23] O. Bergmann, S. Seki and J. Sonnenschein, Quark mass and condensate in HQCD, 
JHEP 0712 (2007) 037,larXiv:0708.2839[ 

[24] A. Dhar and P. Nag, Sakai-Sugimoto model, Tachyon Condensation and Chiral 



symmetry Breaking, JHEP 0801 (2008) 055, arXiv:0708.3233. 

[25] A. Dhar and P. Nag, Tachyon Condensation and Quark Mass in Modified Sakai- 
Sugimoto Model, Phys. Rev. D78 (2008) 06602 1 . larXiv:0804.4807l 

[26] O. Aharony and D. Kutasov, Holographic Duals of Long Open Strings, 
rarXiv:0803.3547! 

[27] K. Hashimoto, T. Hirayama, F. Lin and H. Yee, Quark Mass Deformation of 



Holographic Massless QCD. larXiv:0803.4192i 

[28] K. Hashimoto, T. Sakai and S. Sugimoto, Holographic Baryons: Static Properties 
and Form Factors from Gauge/String Duality, Prog. Theor. Phys. 120 (2008) 1093, 
larXiv:0806.3122l 

[29] M. Edalati, R. G. Leigh and N. Nguyen, Transversely-intersecting D-branes at finite 



temperature and chiral phase transition, arXiv:0803.1277t 
[30] R. McNees, R. C. Myers and A. Sinha, On quark masses in holographic QCD, JHEP 



0811 (2008) 056. iarXiv:0807.5 1271 

[31] P. C. Argyres, M. Edalati, R. G. Leigh and J. F. Vazquez-Poritz, Open Wilson Lines and 
Chiral Condensates in Thermal Holographic OCD. a rXiv:08 11.46171 

[32] D. Bak and H. U. Yee, Separation of spontaneous chiral symmetry breaking 
and confinement via AdS/CFT correspondence, Phys. Rev. D71 (2005) 046003, 
Lhep-th/0412170. 

[33] M. E. Peskin and D. V. Schroeder, An Introduction to quantum field theory, Addison- 
Wesley (1995). 

[34] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products, Fifth 
Edition, Ed. A. Jeffrey, Academic Press (1993). 



25 



